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ABSTRACT 

Weak gravitational lensing by large scale structure affects the number counts of faint 
galaxies through the "magnification bias" and thus affects the measurement of the 
angular two-point correlation function uj(ff). At faint magnitudes the clustering am- 
plitude will decrease differently with limiting magnitude than expected from Limber's 
equation. The amplitude will hit a minimum and then rise with limiting magnitude. 
This behavior occurs because u>(9) due to clustering decreases with distance, while 
the "magnification bias" due to weak lensing increases with distance. The apparent 
magnitude m m i n at which the magnification bias starts to dominate the observed 
clustering is model and color dependent. It is given by u>(m = m m in,9 — 5') f=a 
(1 — 2) x 10~ 3 (5s — 2) 2 Sl§crf , where s is the logarithmic slope of the number counts. 
Already published measurements of to (9) at R = 25 may be strongly influenced by the 
"magnification bias" . An experiment using the ratio of blue and red number counts 
across the sky can be designed such that the effects of the "true" clustering is min- 
imized. The magnification bias is a measurement of the clustering of the mass. This 
weak lensing experiment does not require measuring shapes and position angles of 
galaxies. I derive a revised Limber's Equation including the effects of magnification 
bias. 
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1 INTRODUCTION 

One of the most useful statistics for the evolution and ciustering properties of faint, distant gaiaxies is the two-point angular 
correlation function u>(8). It measures the excess number of pairs of gafaxies separated by angie 6 on the sky. For a given 
galaxy sampfe, lo(6) depends on the redshift distribution of the gaiaxies, and the three dimensional clustering amplitude 
f(r, z). As pointed out by Koo & Szalay (1984), u(6) can be used to constrain models of the evolution and clustering of 
faint galaxies at intermediate and high redshift. A iarge number of investigations have measured the correiation for faint 
magnitudes in different passbands, eg. Bernstein et al. 1994; Brainerd, Smaii, & Moufd 1995 (BSM); Couch, Jurcevlc, & 
Boyfe 1993; Efstathiou et al. 1991; Neuschaefer, Windhorst, & Dressier 1991; Roche et al. 1993. There is a general consensus 
that can be approximated by a power law with slope of —0.8 and an ampfitude that decreases with increasing iimiting 
magnitude. There is, however, iittfe consensus on what that means in terms of fuminosity evoiution, merger history, and 
clustering. The number counts as a function of magnitude, (eg. Liify et al. 1991, Smail et al. 1995), and the redshift surveys 
to faint magnitudes (eg. Broadhurst et al. 1992, Glazebrook et al. 1994) put significant constraints on the galaxy evolution. 

Weak gravitational lensing of distant gaiaxies, equivalent to a systematic induced magnification and distortion without 
multiple imaging, is a probe of the farge scafe mass distribution of the universe, (eg. Gunn 1967). Essentially, when a light 
ray from a distant galaxy passes through an overdense region in mass, the image will be magnified and distorted tangentially 
with respect to the center of the mass overdensity. This effect has been measured in a number of galaxy clusters, [eg. Tyson, 
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Valdes, & Wenk (1990), Bonnet et al. (1994)]. Conversely, passing through a low density region, the image will become dimmer 
and elongated radially. Independent perturbations along the line of sight add up stochastically [Blandford et al. 1991(BSBV), 
Miralda-Escude 1991, Kaiser 1992, Villumsen 1995a(V95a)]. Weak lensing by large scale structure is observable because galaxy 
images near each other on the sky will be sheared coherently leading to a locally preferred direction. Due to the weak nature of 
the lensing there has been no undisputed measure of the distortion, see f.ex. (Mould et al. 1994; Kaiser, Squires, & Broadhurst 
1995), though a reanalysis by Villumsen (1995b) shows a tentative detection at the 5-cr level consistent with expectations 
from a weakly biased flat Cold Dark Matter (CDM) universe. The problem with weak lensing is that it is weak. Measuring 
the weak shear requires measuring the shape and orientation of distant galaxy images. This is observationally feasible, but 
quite difficult, since the images are small and faint, and the possible systematic errors are the limiting factors. 

Weak gravitational lensing is characterized by the convergence k and shear 7. We will work only in the weak limit, i.e. 
k|, I7I <C 1. In this limit, both k, and 7 are observables. The distortion of an image is given by 7, such that an intrinsically 
round image will aquire an eccentricity equal to | -7- 1 and the phase of 7 will be twice the position angle of the major axis. 
The magnification A of an image is A — 1 + 2k. In clusters, the distortion is used to generate shear maps which can then 
be converted into maps of the surface density as pioneered by Kaiser & Squires (1993) and further developed by Seitz & 
Schneider (1995). Recently Broadhurst, Taylor & Peacock(1995) have demonstrated the use of the magnification to improve 
the information about the surface mass density. 

In this paper I propose a new measure of weak gravitational lensing by large scale structure based on the magnification 
bias introduced by Turner (1980) and used by Broadhurst (1995) for lensing in clusters. Basically, weak lensing will induce 
angular correlations in the number density of galaxies on the sky. This measure only involves the convergence n and does not 
involve 7. Thus, it is not necessary to measure the shape of high-redshift galaxies. This means that a galaxy survey to look 
for magnification bias can be pushed to fainter magnitudes than a search for galaxy distortions. 

In §2 I estimate the amplitudes of "true" and "apparent" clustering, in §3 I work out the theory. In §4 I calculate u>(6) 
in terms of the 3-D correlation function and in §5 I discuss the implications. 



2 ESTIMATES OF "TRUE" AND "APPARENT" CLUSTERING 



Magnification bias is the change in surface number density of sources on the sky due to gravitational lensing (eg. Turner et al. 
1984). As the total number of galaxies on the sky is conserved, magnification in a given patch of the sky will lead to a decrease 
in the number density of galaxies, because you are looking at a smaller area of sky than you think. However, magnification 
will also brighten an image so more galaxies will be visible at a given apparent magnitude m. If the number counts No(m) in 
the absence of lensing has a slope 

_ d IogiVo(m) 



dm 

then the number counts will be changed 



(1) 



N oba (m) = No^A 2 - 58 ' 1 = iV (m)(l + (5s - 2)/t), (2) 

where N b s (m) and No(m) are the number counts in the presence/absence of lensing. The first equality is generally valid, the 
second assumes that 7], k| <g; 1. For s = 0.4 there is no magnification bias, for s < 0.4 there is a depletion of counts, while 
for s > 0.4 there is an increase in counts in the presence of magnification. 

In order to estimate the clustering effects of the "magnification bias", assume a universe with no intrinsic clustering of 
galaxies. Then weak lensing will introduce an apparent clustering characterized by an angular two-point correlation function 

([N obB $ + 0,m)-N (m)] [N obs m) - N (m)] ) 
= N^j2 • (3) 

Here the correlation function is calculated as the average over all direction vectors <j> of the relative excess number of galaxy 
pairs separated by angle \6\ = 9. We obtain from Eqs. (^,^) 

( [N (m)(5s - 2Hj + 0)] [N (m)(5s - 2)k{4>)} ) 
UJkk{ ) N (m) 2 { ' 

= (5s-2) 2 (k(4, + 8)k(0)) = (5s-2) 2 C KK (0) = (5s - 2) 2 C PP (8). (5) 

Here, C K k(0), and C pp (8) are the correlation functions of the convergence and the shear. In this weak limit they are identical 
(BSBV and V95a). The shear and convergence correlation functions are steeply increasing functions of the comoving angular 
distance y, and thus of redshift z. For a universe with Qq = 1, J1a = 0, 
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u KK oc C KK oc y s = 8 (l - (1 + z)- 1/2 )\ (6) 

The variation in the number counts to a given magnitude limit across the sky is an observable. The shear and convergence 
correlation functions are related to the density distribution and cosmological parameters as described in BSBV, V95a. The 
shear 7 will also influence the number counts, but the effect is quadratic in 7 and can be ignored. 

In real life, of course, we observe a true clustering of galaxies at faint magnitudes characterized by the angular two-point 
correlation function w M|[x (0), (see references in §1). The correlation induced by magnification bias lj kk (9) is of practical interest 
only if it becomes significant compared to oj mm (#). In a universe where the clustering is constant in comoving coordinates, 
kV M (#) will fall off roughly inversely with distance. Lensing is a cumulative effect while the intrinsic three dimensional clustering 
will be diluted by projection effects. This would indicate that at low redshift u KK (9) is unimportant but that at sufficiently high 
redshift ui KK {9) might become important. The observational issue is whether this occurs at observable apparent magnitudes. 

The current status of ui(6) for red samples is summarized by (BSM), see Fig. 2. It is common to characterize the correlation 
function as a power law with slope 7 and the amplitude at 1 degree, A w . Then lo(9) « A w 9~ 1 , 9 measured in degrees. There 
is a general consensus that 7 ~ 0.8 and that A w is a decreasing function of depth. (BSM) state that the amplitude from red 
counts is A w oc ir°' 27±0 ' 01 with A W (R = 25) « 3 x 1(T 4 . 

For a CDM universe where the characteristic redshift of the sources is z ~ 1, the amplitude of the polarization correlation 
is given by the cosmological density parameter Qo and the rms density fluctuations on a scale of 800 km s _1 , as (BSBV, 
V95a), 

C KK (6 = 5') « 1.5 x 10- 3 Oooi (7) 
This would indicate that the amplitude of A w induced by lensing is 

A w , K &2x 10 _4 (5s - 2) 2 n^|. (8) 

The Mould et al. (1994), presented in (BSM) were obtained in the r band where the source counts have s r = 0.3 . Then 

A w , K (r) « 5 x lO -5 ^ 2 . (9) 

This is lower than the observed value of the correlation function A W (R = 25) « 3 10~ 4 , A W (R — 25.5) « 1.5 10~ 4 , but not 
negligible. However, we have neglected the correlation between the magnification bias and the true clustering. If light traces 
mass, then there will be a strong correlation between the magnification and the true clustering. If s < 0.4, this is really an 
anticorrelation and the observed clustering will be less than the true clustering. 

Smail et al. (1995) find that the slope in the VRI bands all tend to s « 0.3 for the faintest galaxies. Neuschaefer et al. 
(1991) have reported on ui(9) in the Gunn-g band and they find a slope s g w 0.45. This is close to the critical slope of s — 0.4, 
so the magnification bias is expected to be small and quite uncertain. For g < 24.5 they find that lj(9) is approximately a 
power law with the amplitude decreasing with increasing magnitude. For fainter magnitudes they find a very steep rise in the 
amplitude. This rise is much too steep to be accounted for by magnification bias. At g w 24.5, ui(6 = 1') w 10~ . For a high 
value of f2o crs there may be a small contribution to lj(9) from magnification bias, though this is uncertain. 

However, Broadhurst (1995) has shown that for the / band counts, / > 24, the red counts, V — I > 2.0, are fitted well 
by a slope s r w 0.15. The blue counts, V — I < 1.0, have s b w 0.5 . The prediction is thus that for the blue counts, there will 
be a measurable magnification bias. For the red counts, though, the prediction is that the magnification bias will be strong, 

A w , K (r) w 3 x 10- 4 n 2 al. (10) 

This is comparable to the observed correlation function amplitude. 

These results indicate that at accessible apparent magnitudes the magnification bias will give a significant contribution 
to the observed angular two-point correlation function. This encourages a more detailed study of the magnification bias. 



3 THEORY 

3.1 Single Sample Statistics 

In the limit where the intrinsic clustering and the magnification bias are weak, the clustering, characterized by the angular 
two-point correlation function u)(6) can be calculated. The magnification bias works on the "true" number counts, i.e. the 
number counts which include the true clustering, 

N oba (m) = JV (m)(l + A^A 2 53 ' 1 « N (m)(l + A/i + (5s - 2)k). (11) 

Here, A/i is the true fractional excess number of sources in a particular patch on the sky. From this we can calculate ui(9) 
from an average over the sky, 

( [N obs $ + 9, rn) - No (m)] [N obs m) - N (m)] ) 
iVo(m) 2 



u(fi) = m ^1 : >_ — , — (12) 
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= ( [A/i(<£ + 9) + (5s - 2)k{4> + §)] [A/i(<£) + (5s - 2)k(0)] ) (13) 
= (Aii(4 + 8)Aii($)) + (10s - 4) (An{$ + 3)k$)) + (5s-2) 2 ( k (0 + 0) k (0)) (14) 
= w MM (6») + 2w MK (0)+o; KK (6l) (15) 

= u J ^(6) + 2u JllK (e) + (5s~2) 2 c KK (e). (16) 

Again, the average is over all direction vectors (j>. The three individual terms are not observables, only the sum. The first 
term is the true galaxy clustering, the third term is the apparent clustering induced by the magnification bias, and the second 
term is a cross correlation. The magnification bias is induced by mass density fluctuation, while the true galaxy clustering 
measures the galaxy number density fluctuations. If in projection, the number density fluctuations are decoupled from the 
mass density fluctuation, then the cross term will be zero. 

The observed number counts of galaxies as a function of magnitude in a given passband and position on the sky will 
depend in a complicated, and unknown way, on the formation and evolution of galaxies. In the calculations we are going to 
use the comoving radial distance x as fundamental variable, 

x(z) = f dz' H- 1 (z') ; H(z') = [^(1 + z'f + (1 - fi„ - Q A )(1 + z') 2 + fi A ] 1/2 , (17) 
Jo 



sinh 

y{z) = 



- n - n A ) x{z) 



(18) 



In a simple model, characterised by a selection function S(x), and linear bias factor b(x) we can calculate statistically the 
surface density distribution. We define S(x) to be the mean comoving number density of observed objects at distance x. The 
function S can also be thought of as a function of redshift, or lookback time. The normalisation of S is such that 



f 

Jo 



dx y 2 S(x) = 1. (19) 

/o 

The function S has hidden in it the number density evolution, the luminosity evolution, and the differential K-correction. We 
define the linear bias factor b as 

*>-(#)/(?)■ 

so b is not a function of scale but can be a function of epoch. This linear biasing scheme assumes that the galaxies cluster the 
same way as the mass, just at a possibly different amplitude characterized by b. 

The surface density fluctuation A/li can then be written as a line integral in comoving radial distance x of the fractional 
volume overdensity b(x)8(x) times the volume element dx ■ y 2 with a selection function S(x). Here y is the comoving angular 
diameter distance. 

POO 

A/j(0) = / dxy 2 S(x)b(x)8(x), (21) 
Jo 

/•oo 

k{9) = 3O / dx yw{x)5(x) / a. (22) 
Jo 

In principle, all the line integrals are only out the horizon distance xh, however, we can just set S(x) equal to zero for x > xh- 
The equation for k is taken from V95a, Eq 29, where n = —AM. Here, 5(x) is the density perturbation at position x evaluated 
at the epoch when the lightray passes through. The expansion factor a also enters, a = 1 at the present epoch. The function 
w(x) is the lensing selection function which for a lens at distance x is the integral over all sources y 2 S(x') further away than 
x of the ratio of the lens-source distance yLS and the observer-source distance yos- 

p oo 

w(x) = / dxy 2 S{x)^- (23) 

For convenience we write 8(x) = f(x)aSo(x), where 5o(x) is the density fluctuation evaluated today. We implicitly assume 
that the "growth factor" f(x) is a universal function of expansion factor a(x) only. Then 

/•oo 

A/x(0) + (5s - 2)k(4>) = / dxy f(x)[y S(x)ab(x)+3n (5s-2)w(x)]5 (x). (24) 
Jo 

Suppose we have sufficiently well behaved functions F(<j>), and G(x) defined as 

f X H 

F(4>) = / dxG{x)5 (x). (25) 
Jo 

The two-point correlation function Cff{0) of F will then be 
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C FF (0)=4n dxG 2 (x) dkk P(k)J (kx9). (26) 

Jo Jo 

Equation [m| is of this form so we can write the two-point correlation function as 

uj{9) = 4tt 2 / dx y 2 f 2 (x) [yS(x)a b(x) + 3Q (5s ~ 2)w{x)] 2 / dkkP(k)J (kx9). (27) 



We see in Eq.(^) that if s > 0.4, the weak lensing will give a positive contribution to the correlation amplitude, while if 
s < 0.4, the weak lensing will have a negative contribution to the correlation amplitude. 

The weak lensing can change the angular dependence of uj{9) if P(k) is not a power law. On a given angular scale 9, to(9) 
samples decreasing wavenumbers as a function of x as seen in Eq.|27]. The weighting in k is given by the outer integral which 
is changed by the magnification bias. Thus, if P(k) is not a power law u) (9) will be changed by the magnification bias. 

At low redshift the intrinsic clustering will dominate and the magnification bias is negligible. At higher redshift, the 
intrinsic clustering will drop and the magnification bias will increase dramatically. In the case of flat number counts, i.e. 
s < 0.4, the angular two-point correlation function will decrease faster with depth than expected from Limber's equation 
(Peebles 1980). Eventually the contibutions from the true clustering and the magnification will be about equal and the 
observed clustering amplitude will hit a minimum. At even higher redshift the magnification bias will dominate and uj(9) will 
increase. 

For steep number counts, i.e. s > 0.4, the correlation amplitude will fall slower with apparent magnitude than expected 
from Limber's equation, reach a minimum, and then rise again. In both cases the correlation amplitude will eventually become 
an increasing function of depth. The magnitude at which this happens will depend on cosmological parameters, Qq and b. 

How do the intrinsic correlations and the apparent correlations scale with depth of the sample. Take a simple example 
where S~lo = 1, =0, 1 - a < 1, and b{x) is a constant, then y = x. In Eq.(|2^) in the outer integral, the contribution from 
intrinsic clustering will scale inversely with q, while the contribution from magnification bias will scale as q A . In Eq. ^ the 
calculation of the intrinsic clustering involves five powers of the distance. As S 2 is inversely proportional to the sixth power 
of the characteristic distance of the sources we have that ^^^(9) scales inversely with characteristic distance of the sources. 
Also in Eq. |27] we see that ui KK (9) scales as the third power of the characteristic distance. The inner integral in k is the same 
for the two terms, and is proportional to x~ 2 ~ n for a power law power spectrum with slope n. 

We can thus estimate the sample distance at which the magnification bias will dominate. The correlations induced by 
weak lensing are approximately uj KK {9) oc f^o-fy 1- " while the true angular correlation has 10^^(9) tx a 2 b 2 y~ 3 ~ n . Thus the 
characteristic distance y e at which the two effects become comparable is y^ 2 oc Qo/b. The limiting magnitude m m i„ at which 
the correlation amplitude is at a minimum is an observable. If the typical sample redshift is z ~ 1 then m m i n is determined 
from the equation. 

u(m = m min ,6» = 5') « (1 - 2) x 10~ 3 (5s - 2) 2 fio<7 2 . (28) 

In Figure 1 we demonstrate the effects of the magnification bias through a simple model. Assume that S(x) is constant out 
to some distance xq. Assume that b = 1/og, i.e. the bias is independent of epoch. Further assume Qo = 1 and P(k) oc k~ 12 . 
In this model, which should be seen as an illustration only, we can calculate Aw as a function of mean source redshift 
< z >= f Q dx z(x)x 2 S(x), for various values of s and b. In all these models the intrinsic clustering is the same. 

The solid curve, s = 0.4, is due entirely to the intrinsic clustering, there is no magnification bias. The uppermost curve, 
s — 0.5, b — 1, shows a positive magnification bias and falls off slower than the true clustering. At < z >~ 1 the observed 
clustering amplitude is 70% higher than the true amplitude. At higher redshifts A w is quite flat and at < z >~ 2.5 it will 
begin to rise. 

The two curves for s = 0.3, which are appropriate for red samples, originally fall faster than the true correlation. For 
z Si 1.8, the less biased model, dotted line, will have a lower amplitude of A w than the more biased model, while at higher 
redshifts the opposite is true. This behavior can be understood as follows. The magnification bias consists of two terms, u KK (9) 
which is always a positive contribution and quadratic in <rg, and w (JK (0) which is linear in erg and is negative. At low redshift 
the negative linear term dominates while at sufficiently high redshift the positive quadratic term dominates. 

For an unbiased population of galaxies with s = 0.2, the magnification bias is so strong that the correlation amplitude 
hits a minimum at < z >~ 1.3 and is larger than the true amplitude for < z > <; 1.5. The curve for b — 2, s = 0.2 is identical 
to the curve for b = 1, s = 0.3. 

If we observe the correlation function u(9) for a single limiting magnitude we cannot tell what is the intrinsic contribution 
and what is magnification bias, it is necessary to use a range of limiting magnitudes. The slope of the number counts s is 
an observable, so if s < 0.4, and the correlation amplitude is decreasing with limiting magnitude, we know that the intrinsic 
correlation is higher than the observed correlation. If the amplitude is rising we cannot make this inference. If s > 0.4 we 
know that the true correlation is less than the observed correlation. 
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Figure 1. Correlation Amplitude A w as a function of mean redshift and slope of number counts. The solid curve is the intrinsic 
correlation. The dotted, long dashes, dash-dotted curves are relevant for red galaxy samples. The short-dash curve is relevant for blue 
galaxy samples. 



If an increase in correlation amplitude as a function of limiting magnitude is observed, we see that it is not necessary to 
invoke a local population of weakly clustered intrinsically faint galaxies. It is also not necessary to assume that we are seeing 
a strongly clustered population of distant galaxies. 

These models should be seen as an illustration only. In order to properly model the correlation function we need to 
consider luminosity evolution, merger history, density evolution, differential K-correction etc. In particular we need to relate 
< z > to limiting magnitude. 

From Figure 1, we see that if the data at R w 25.0 — 25.5 in (BSM) is at < z >~ 1, there is likely to be a significant 
contribution from magnification bias in the data. 



3.2 Two-Sample Statistics 

An equivalent statistic involving two samples can be used. The number counts of sources selected in two separate ways in the 
sample field are used. The ratio of counts in two passbands at magnitudes (mi, 7712) can be calculated in the case of weak 
intrinsic clustering of galaxies. 

A , 7, N obs (mi,(j}) / N obs (m 2 ,<j>) 1 + Ap(mi, $) + [5si - 2]ki(0) 

[ mi ' 2,<W " JVo(mi) / N (m 2 ) ~ 1 + A M (m 2 , $) + [5s 2 - 2]« 2 (0) 1 ' 

w A^(mi,0) - Afi(m 2 ,<j>) +5 [si/ti(0) - s 2 n 2 (4>)] - 2 - K 2 (<j>)] . (30) 
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Figure 2. Correlation Amplitude A w as a function of mean redshift and difference in slope of number counts. The three curves that 
start below the solid curve arc for s\ < S2- The three other curves are for s\ > S2- 



(7711,7712), (si,S2), (ki,K2) are the magnitudes, slopes of the number counts and the convergences in the two passbands. The 
convergences need not be identical since the depth of the samples in redshift need not be the same in the two passbands. 
Equivalently, we can look at two separate samples selected by color, in a single passband. This statistic in the weak clustering 
regime is equivalent to the statistic involving the difference in number counts. 

N h s (mx,4>) N oba (rri2,cj>) 



A(77U, 7712 



TVo(mi) No(m 2 ) 
= A/i(m 1 ,^) - Afi(m 2 , <j>) + 5 [si«i(</>) - S2K2WO] - 2 [«i(0) - k 2 (<^)] 
The correlation function can then be calculated in a similar fashion to Eq. ( pTj ) 



(31) 
(32) 



u>{9) 



4tt 



dx y f (x) [y a (Si(x)bi(x) — S2(x)b2(x)) + 3ilo (5(sitvi(x) — S2W2(x)) — 2 (wi(x) — W2(x)))Y 



dk k P(k)J {kx6). 



(33) 



Figure 2 shows the correlation amplitude for the same model model as in Figure 1. Here we have assumed two populations 
with the same selection function Si(x) = S2(x), thus wi(x) = W2(x). We have further assumed that as = 1 and that 
61 = 1, 62 = 1/2. This means we have strongly and weakly correlated populations. We show A w for As = si — S2 = 
(—0.3,-0.2,-0.1,0,0.1,0.2,0.3). The solid curve is A w in the absence of magnification bias As = 0. For a positive value of 
As the curves are always above the solid curve. At low redshift there is little magnification bias and the curves all follow the 
solid curve. For positive As the curves move significantly away from the solid line at < z >~ 1/2 and then become nearly flat 
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out to < z >~ 2. For negative As, the curves trace the solid line more closely, but with a significant difference at < z >~ 1/2. 
For the curves with a substantial difference in number counts slope, the curves then hit a minimum and then rise steeply. 

These curves for negative As may be relevant for a strongly clustered red population with small number counts slope 
and a weakly clustered blue population with a larger value of s. 

This estimator is most useful when the intrinsic clustering is similar in the two bands, while the slopes of the number counts 
are very different. In the case where the intrinsic clustering is the same in the two passband, i.e. A/i(mi,0) = A/j,(m 2 ,<f)), 
there is zero contribution from the intrinsic clustering. If also the depth in comoving distance is similar, then ki « k 2 and we 
get that 



uw(0) = 25(s 6 - s r ) 2 C KK (9). 



(34) 



This assumption is not as severe as it might seem as can be seen from the following example. Look at two source planes at 
distances x 1 < x 2 and two lightrays separated by a small angle 9. A mass fluctuation at distance x will coherently influence 
the convergence in the two beams only if the beams both pass through the density fluctuation. This is obviously only possible 
if x < x±. Thus density fluctuations beyond the nearest source plane will not influence the correlation, and 



+ 6)k 2 (§)) « (ki(4> + 9)ki[9)) => C K1K2 (6) » C K1K1 (6). 



(35) 



In that case we have to a good approximation that the observed two point correlation function of the ratio of counts is given 
by Eq. [$4| which for the ratio of very red and very blue samples gives 



uj kk (9)^3C kk (9). 



(36) 



If at intermediate magnitudes we find that the observed correlation strengths differ by a factor C 2 we can use a different 
estimator of the clustering, 



A(mi,ni2, 



N obs (mi , 4>) 



C 



N oba (m2, < 



(37) 
(38) 



JVo(mi) N (m 2 ) 
= An(mi,$) - CA/i(m 2 ,<£) + 5 [siKi($) - C s 2 k 2 ($)] - 2 [m(<t>) - C k 2 {4>)] 
The correlation function can then be calculated in a similar fashion to Eq. ( p^] ) 

/>oo 

uj{6) = 4tt 2 / d^y 2 / 2 (x)[ya(S , i(x)6i(i)-CS , 2 (a;)&2(x)) + 3a)(5(simi(x)-C , S2w;2(a;))-2(mi(x)-C*TO2(x)))] 2 
Jo 

x / dkk P(k)J (kx9). (39) 
Jo 

By construction, the intrinsic correlations cancel out at intermediate redshift. At higher redshifts this cancellation may not 
occur, but this estimator may minimize the effects of the intrinsic clustering even at faint magnitudes. 



4 CORRELATION FUNCTIONS, REVISED LIMBER'S EQUATION 

Suppose we have sufficiently well behaved functions F\{<i>), F 2 (4>), 5i(x), 5 2 (x), Gi(x), and G 2 (x) defined as 



dxGi(x)6i(x) ; F 2 {<j>) 



dxG 2 (x)6 2 (x) ; F{4>) = F 1 {4>) + F 2 (<j>). 



(40) 



The two-point correlation function Cff(9) of F will then be 



C FF {9) = 2 



dx 



G\{x) dq Z 11 ((M) 2 + q 2 ) 



2\l/2 



+ 



Gi(x) dqt 



((x9f+q 2 ) 



2\l/2' 



2G!(X)G 2 (X) / dql; 



((X9f + q 2 ) 



2\l/2 



(41) 



where f 11 , £ 22 , and ^ 12 are the two autocorrelation functions of 5i(x) and 8 2 (x) and their cross correlation function. This 
follows from Eq. || and V95a, Eqs. 37-39. 

Let us take some galaxy density distribution S 9 (x) and mass density distribution 5 m (x). Then the surface density of 
galaxies A/i and convergence k are 



Afi(9) = 



dx y S(x)S 9 (x), 



(42) 
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/>oc 

k(6) = 3f2o / dx yw(x)8 m (x) /a. 
Jo 

These equations are useful for calculating ui(6). We obtain 



(43) 



dx 



y 4 S 2 (x) / dg£ 



(M) 2 + 



2\l/2 



+ 



9n () (5s — 2) j/ iu (x)a 



dqC 



((x6) 2 + 



2\l/2 



+ 



6tt (5s-2)y 3 S{x)w{x)a~ 1 dq £ 9m ((:r(9) 2 + g 2 ) 



2\l/2 



(44) 



Here £ mm , and £ 9 *™ are the galaxy-galaxy, mass- mass, and galaxy-mass correlations evaluated at the appropriate epoch. 
In the case where s = 0.4, i.e. no magnification bias, this reduces to the ordinary Limber's equation. 

We can rephrase the correlation function calculation as an integral in redshift. This becomes simply 



lo(6) = 2 1 dz 



H(z)N 2 (z 



dq e 



((xdf + q 



2\l/2 



//^(z) 9fig(5s - 2) 2 y 2 w 2 (z) (1 + zf 



6fi (5s-2) y N(z)w(z) (I + z) / dq £ 



+ 



((x6) 2 + < 



+ 



((x9) 2 + 



2\l/2 



(45) 



Here N(z) — y 2 S(x)H^ 1 (z), w(z) = w(x). Notice that the normalisation integral for N(z) still applies, (Eq. |is| ). From this 
equation it is possible to calculate to(0) given a cosmological model i.e. Qo, ^a, and a selection function N(z) if you know, or 
assume, the three correlation functions. 

Let us make life easy and assume that the correlation functions are all power laws, not necessarily with the same slopes 
and evolution. For each of them assume that in proper coordinates r 

-7 



£(r,*) = 



(1 + 2) 



-(3+«0 



(46) 



Here, the exponents 7 and e need not be the same for the correlation functions. Even more importantly, the proper correlation 
lengths ro need not be the same. With this assumption 



u(0) = lu 39 {6) + Lu mm (9) + 2 Lo 9m (8). 



L0 39 (6) = 



\e) = 



,1 — 1 



uj 9m {6) 



r[( 7 -i)/2] 
r[7/2] 

r[(7-i)/2] 



r[7/2] 



3 f2 (5s - 2) rl 



dz H(z) N 2 {z) x 1 -" 1 (1 + z) 7-3-E , 

dz H~\z) w 2 (z) y 2 x 1 - 1 (1 + z) 7_1_E , 
dz N{z) w{z) y x 1 " 1 (1 + Z y- 2 -\ 



7 zji — r 



(47) 
(48) 
(49) 
(50) 



The term for the intrinsic clustering u gB reduces to the same result as (BSM), Eqs. (6,7) when taking into account 
differences in notation. Let us make life even easier for ourselves and assume that 7 is the same for all three correlation 
functions. However allow e to be different so that e m,n = e + Ae, and e gm — e + Ae/2. Further assume a linear bias model so 
that today = 0) = b 2 £ mm (z = 0), and £ 9m (z = 0) = b ^ 9m (z = 0). Now r is the correlation length of the galaxies. This 

simplifies the integral for u>{&). 



u(6) 



r[(7-i)/2] 
r[7/2] 

dz H(z) 



.7 a 1 — 1 



N(z) + 3 ^ (5s - 2) w(z) y(l + z) 1 '^ 2 H'^z) 



1-7 



(1 + *) 



7— 3— e 



(51) 



Parametrized this way it is possible to do an analysis of uo{6) in the same way as BSM while including the effects of the 
magnification bias. For a given magnitude limit, assume a cosmological model, i.e. Qo, ^a- Then choose 7, e, Ae, b, and ro. 
For a given N(z) this will then predict the observed uo(9). 

As stated in BSM, for 7 = 1.8, linear theory predicts that e = 0.8, clustering fixed in comoving coordinates will give 
e = — 1.2, while clustering fixed in proper coordinates will have e = 0. 

If we have further information about n we can improve on our measurements. If we can measure image shapes and position 
angles, then we can infer the gravitational shear field p. We can then measure the two point correlation function C pp (8). In 
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the weak limit, C KK (6) = C pp (0) and u mm (0) = (5s - 2) 2 C KK (8). Thus we can remove uj mm from our measurements. This 
leaves us with the true clustering of the galaxies plus the cross term galaxy-mass. If the observational data is of sufficiently 
high quality, we can estimate p at any given point on the sky. This estimate is of course smoothed over a finite area. If 
we can measure p across the sky, we can infer k through an inversion procedure such as demonstrated by Kaiser & Squires 
(1993), or Seitz & Schneider (1995). With these methods we can obtain an unbiased estimate of k. There is no problem with 
non-linearities since we are certainly in the weak limit and there is no ambiguity with the mean surface density since it is by 
assumption zero. If we know k then we simply subtract the term (5s — 2)n(cj>) x No(m) in Eq. ^. We thus retain an unbiased 
estimate of the true galaxy clustering. 

With this measurement we can thus separate the intrinsic clustering of galaxies, the clustering of the mass, and the 
galaxy-mass correlation. In other words we can separate the galaxy clustering evolution from the mass clustering evolution, 
and we can also measure how well light traces mass. 

5 DISCUSSION 

The standard way of measuring the galaxy number density fluctuations at intermediate redshifts z ^ 1 is through the angular 
two-point correlation function uj{0). Measuring U)(9) is in principle straightforward. You measure the positions and magnitudes 
of the galaxies and measure the number of galaxy pairs relative to random. Small distortions in the telescope optics do not 
influence the results and it is not necessary to measure the shape and orientation of the galaxies. The interpretation of u>(8) 
in terms of the three-dimensional correlation function £(r) is complicated through largely unknown luminosity and density 
evolution. Even if £(r) were known accurately, this would only tell us the galaxy distribution, not necessarily the mass 
distribution. The mass distribution would then have to be inferred through a model dependent biasing scheme, or some other 
modeling scheme. 

Weak lensing by large scale structure is a direct measure of the mass distribution and it circumvents inferring the mass 
distribution from the galaxy distribution. The problem with weak lensing is that it is weak. It requires measuring the shape 
and orientation of faint galaxy images. This is observationally feasible, but quite difficult, since the images are small and faint, 
and the possible systematic errors are the limiting factors. However, the payoff is immense in terms of measuring cosmological 
parameters and the power spectrum of density fluctuations. 

We have presented a new way of measuring the mass density fluctuations at intermediate redshifts by measuring the 
angular two-point correlation function u>(6). At sufficiently faint magnitudes, the observed u>(8) will be dominated by weak 
lensing. The method combines the relative ease at which u(9) can be determined, with the relatively simple theoretical 
interpretation of weak lensing in terms of cosmological parameters and the statistical properties of the mass distribution. 

At relatively bright magnitudes, R 5s 23, weak lensing is expected to be unimportant. If the slope of the number counts is 
greater/less than 0.4, the amplitude of lo{6) will decrease slower/faster than expected from Limber's equation. In both cases 
the amplitude will hit a minimum and then increase with limiting magnitude. 

An equivalent measure is the ratio of number counts in two different passbands as a function of position in the sky. This 
measure can, by carefully choosing the passbands and magnitudes, be tuned to be nearly independent of the true clustering 
of galaxies. In that case, u){6) is a straight measure of weak lensing and the theoretical interpretation is considerably simpler. 

In summary: The angular two-point correlation function of galaxies is affected by weak lensing by large scale structure 
through the magnification bias. At faint magnitudes this can be a signicant effect and must be included in calculations. For 
a properly designed experiment, ui(9) can be used to infer the clustering of the mass, the clustering of the galaxies, and how 
well light traces mass. 

More work needs to be done. 
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